Depinning of semifiexible polymers in (1 + 1) dimensions 
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We present a theoretical analysis of a simple model of the depinning of an anchored semifiexible 
polymer from a fixed planar substrate in (1 + 1) dimensions. We consider a polymer with a discrete 
sequence of pinning sites along its contour. Using the scaling properties of the conformational 
distribution function in the stiff limit and applying the necklace model of phase transitions in 
quasi-one-dimensional systems, we obtain a melting criterion in terms of the persistence length, the 
spacing between pinning sites, a microscopic effective length which characterizes a bond, and the 
bond energy. The limitations of this and other similar approaches are also discussed. In the case 
of force-induced unbinding, it is shown that the bending rigidity favors the unbinding through a 
"lever-arm effect" . 
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I. INTRODUCTION 

For a broad category of physical problems, a free poly- 
mer is characterized by two lengths: the total contour 
length (L) and the persistence length (L p ) which is the 
correlation length of the tangent unit vector along its con- 
tour and is proportional to its bending rigidity. When the 
persistence length is much smaller than the total length, 
the polymer is said to be flexible and can be treated as 
a random walk. When the two lengths are of the same 
order, the polymer is said to be semifiexible. Some of the 
most important biopolymers belong to the latter class. 
For example, the structural elements of the cytoskele- 
ton are microtubules, actin filaments, and intermediate 
filaments with persistence lengths of the order of 6mm 
0, 17^m || , and 2^m g] respectively. Although DNA 
filaments usually have a total length greater than the per- 
sistence length (L p w 50nm), the latter is long enough to 
affect their elastic properties. (J] Obvious biological rel- 
evance and inherent theoretical challenges have sparked 
great interest in the statistical mechanics of semifiexible 
polymers in recent years. Q 

A theoretical analysis of the unbinding of semifiexible 
polymers from fixed surfaces or interfaces (adsorption- 
desorption transition) or of two semifiexible strands from 
each other is a particularly tricky problem. The main 
reason is that sharp phase transitions in statistical me- 
chanics occur only in the thermodynamic limit and the 
thermodynamic limit of semifiexible polymers is ambigu- 
ous. If we keep the persistence length fixed and take the 
total contour length to infinity, we obtain a flexible poly- 
mer. If we take the persistence length to infinity keeping 
the total length fixed, we obtain a rigid rod without any 
fluctuations. There have been several studies of this sub- 
ject over the past few years. § g, §, g, |T| [l2| [D§ In 
all of those works the polymer binds to a potential well 
which continuously extends over the surface (or inter- 
face). In this article, we consider a simple model where a 
weakly bending semifiexible polymer in (1 + 1) dimensions 
is bound to a fluid surface (or interface) through a dis- 



crete sequence of regularly placed pinning sites ("sticky 
points") along its length. Such a model resembles the 
physical situation where an actin filament binds to a 
membrane through anchoring proteins. fl4|| The discrete- 
ness of the binding sites allows us to employ a different 
kind of thermodynamic limit which avoids the inconsis- 
tencies that appear in some previous works. For fixed 
total length and persistence length (with L <C L p ), we 
take the density of binding sites to infinity under the con- 
straint that the probability of finding such a site inside 
the binding region remains constant. Thus we obtain a 
second order unbinding transition and a melting temper- 
ature which is a function of of the persistence length, 
the spacing between pinning sites, a microscopic length 
which characterizes a bond, and the bond energy. The 
article is organized as follows: In Sec. II we calculate 
the probability of finding a binding site of the polymer 
inside a small (microscopic) region which characterizes a 
bond. We then consider the simplest version of our model 
which is a filament with only one pinning site. In Sec. 

III we use a necklace model |l6| type of approach 
which yields the thermal depinning transition. In Sec. 

IV we discuss the effect of the bending rigidity on the 
force-induced unbinding of semifiexible polymers in the 
stiff limit. Finally, in Sec. V we demonstrate the sub- 
tleties of this problem comparing our model with other 
approaches and we present our conclusions. 



II. CONFORMATIONAL PROBABILITY - 
FORMALISM 

A widely used model which captures much of the 
physics of semifiexible polymers (except for their self- 
avoidance) is the worm-like chain (WLC)|l7]] where the 
polymer is considered to be a continuous inextensible 
curve r(s) parametrized by the arc length s measured 
along its contour from a fixed end. The effective free en- 
ergy of a particular conformation depends only on the 
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bending (curvature) and is given by 
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where t(s) = dr(s)/ds is the tangent unit vector of the 
curve r(s) and n is the bending rigidity which is related 
to the persistence length via L p = 2n/ksT (in two di- 
mensions) . 

The orientational probability distribution function for 
free semiflexible chains having an initial tangent vector 
t(0) = to and a final tangent vector t(L) = is given 
by the path integral 



G(ti,L|t ,0) = 

r t(L)=t L 

Af / V[t(s)}S[\t(s)\ 
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where the integration is over all fluctuating "paths" t(s) 
subject to the fixed boundary conditions and the inex- 
tensibility constraint |t(s)| = f. TV is a normalization 
constant. There is a formal analogy between the clas- 
sical statistical mechanics of a semiflexible polymer and 
the quantum statistical mechanics of a rigid rotator. [p~7| 
If we make the correspondence k ^ /, fc^T ^ h, and 
L ^ f3h in Eqs. (g,|), we notice that G(t L ,L|t ,0) 
corresponds to the density matrix element, in the angle 
representation, of a quantum rigid rotator with moment 
of inertia I and inverse temperature (3. As in the case 
of a density matrix, ]l8|] the angular probability distri- 
bution function of a free semiflexible polymer satisfies a 
Schrodinger equation in imaginary time: 



dG _ 1 d 2 G 
ds ~ L p d9 2 



(3) 



where 9{s) is the angle between t(s) and a fixed refer- 
ence axis. 19 In order to obtain the complete distribu- 
tion function which in addition to the tangent vector also 
includes the position vector, G(r s ,t s , s|r , t , 0), we have 
to replace the s-derivative in the lhs of Eq. (||) by the 
"convective" derivative d s + t ■ V r along the polymer 
"path" r(s) with instantaneous position vector r and tan- 
gent vector t.[^o) In Cartesian coordinates, the equation 
reads 
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G(x a ,y s ,e s ,s\x ,yo,e ,0) =0 (4) 



where 9 is the local slope of the polymer with respect to 
the x-axis. 

In the weakly bending limit (L <C L p ), 9 -C 1 and we 



simplify Eq. (|J) setting sin# » 9 and cos 9 ~ 1. Since 
we are not interested in the longitudinal fluctuations of 
the polymer (along the x-axis), we integrate the com- 
plete probability distribution function over x to obtain a 



simpler equation for the reduced probability distribution: 
d 1 d 2 



d_ 
ds~ 
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G(y s ,9 s ,s\y ,9 0l 0) = . (5) 



Using Fourier transformations |21| 
the "initial" condition lim s .^o G(y s 
0o)8(y - yo) to get 

G{y a ,e s ,s\y ,e ,0) = 



we solve Eq. ( 
Sa,s\y ,9 ,0) 
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Apart from explicitly containing the persistence length 
L p , Eq. (H) is identical with that obtained in Ref. [0. 
The interpretation, however, is very different. In Refs. 
H 0], G(y s , 9 S , s\yo, 9q, 0) is interpreted as a dimen- 
sionless partition function independent of the persistence 
length L p which has been eliminated by rescaling y and 
9. In those references, Eq. ^ is expected to be valid 
for large s and it appears that s is measured in units of 
an extra, "monomer" length. In contrast, we interpret 
it as a two-point conformational probability distribution 
valid only in the weakly bending limit (s -C L p ). No- 
tice that G(y s , 9 S , s\yo, 9q, 0) fulfils the three fundamen- 
tal properties of a two-point probability distribution; its 
integral over y s and 9 S is 1, it becomes a delta func- 
tion when s — > and it obeys the Chapman-Kolmogorov 
equation. p6| The corresponding partition function dif- 
fers from G(y s , 9 S , s\yo, 9q, 0) by a normalization fuctor 
(related to the measure of the path integral) which should 
have units of length in order to render it dimension- 
less. (It is similar to the phase volume element 2ir?i 
used in the statistical mechanics of gases.) In the cal- 
culation of several quantities, this normalization factor 
is unimportant as it drops out. For this type of prob- 
lems, G(y s ,9 s , s\yo,9o,0) itself can be considered as the 
partition function. However, as it will become clear be- 
low, the necklace model involves a sum over powers of 
the partition function and using a dimensionful quantity 
in its place would clearly be erroneous. 

For fixed yo = 9q = 0, the mean square slope and 
transverse displacement of the free end of a filament of 

as 



length L are (ffj) = 2L/L p and (yV) = (2/3)L 3 /X : 
can be easily calculated from Eq. (o). The probability of 
finding the free end within a very small range of slopes 
and transverse displacements (—6 < 9l < 5 and — e < 
y L < e with < 5, e < 1) is 

P{S,e,L,L p ) = J dO L J dy L G{y L ,9 L ,L\0, 0,0) 



2^ B 



(7) 



where B = 4<5e and the approximation holds for B <C 
(V3/2tt)L 2 /L p . The partition function Z(S, e, L, L p ) of 
a polymer which is constrained so that yo = 9q = and 



— 8 < 8l < S , — e < yi < e while it is unconstrained 
in the longitudinal direction is related to the probability 
P(<5, e, L, L p ) via Z{5, e, L, L p ) = Zf{L, L p )P{5, e, L, L p ), 
where Zf{L,L p ) is the partition function of a free fila- 
ment. The latter has the property Zf(L\ + L 2 , L p ) — 
Zf(Li,L p )Zf(L2,L p ) and will be neglected as it is not 
going to affect any of the observable quantities we are 
interested in. 

The probability of finding both the free end and the 
point at the middle confined within a very small range of 
slopes and transverse displacements is 

/ Ml J dy L G(y L ,e L ,L\y L/2 ,e L/2 ,L/2)x 
s 
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We shall use this factorization in the calculation of the 
partition function of our model. 

We now consider the toy system of a weakly bend- 
ing semiflexible polymer with its endpoints (s = and 
s = L) bound and a pinning site in the middle (s = L/2). 
A "bound site" in our model is defined as a point of the 
polymer which is constrained to fluctuate within a micro- 
scopically small range of slopes and transverse displace- 
ments which is characterized by the effective length B as 
defined above but is free to fluctuate in the longitudinal 
direction. The latter situation is physically realized in 
the case of a fluid substrate where the "sticky points" 
are free to move along a Id track (membrane). A "pin- 
ning site" is defined as a point on the polymer which is 
energetically favorable to be bound with an associated 
bond energy J (J > 0). The partition function of this 
system with one pinning site is 



where 



Z 1 = [g{L/2)] 2 v + Q{L)-[Q{L/2)f 



V3 L p 



(9) 



(10) 



is the conformational statistical weight of a polymer seg- 
ment of contour length I whose endpoints are bound and 
V = exp(J/k B T) . The third term in Eq. (|) is the 
"counterterm" needed to prevent doublecounting of con- 
formations; the conformations associated with [Q(L/2)] 2 
have already been included in Q{L). 
The average fraction of intact bonds is 



Q = 



dlnZi 
<91n(v - 1) 



(11) 



This is a general expression valid for any number of pin- 
ning sites provided that we replace Z\ with the corre- 
sponding Zjq and we divide the rhs by N . The calcula- 
tion of Zm for N ^> 1 is the aim of Sec. III. 




FIG. 1: A cartoon picture of a polymer conformation with a 
"chain" of length 3L m , a "bubble" of length 6L, n , a "chain" of 
length 2L m and a "bubble" of length 6L m . The dots represent 
pinning sites and the vertical black lines represent anchoring 
proteins which bind the filament to the substrate. The weakly 
bending limit allows us to neglect any direct (hard wall) in- 
teraction with the substrate. 



III. THERMAL DEPINNING TRANSITION 

The partition function of a weakly bending semiflexible 
polymer with its endpoints (s = and s = L) bound 
and N pinning sites regularly distributed along its length 
formally reads 



Zn = 



N 

E 

n=0 



v n V(n) , 



(12) 



where V(n) is the probability of a conformation with ex- 
actly n bonds {but not n + 1 or n + 2 or... N bonds). 
For example, in the case of N = 2, V{2) = [£(L/3)] 3 , 
V{1) = 2£(L/3){S(2L/3) - [G{L/3)} 2 }, and V{0) = 
G{L) - 2G{L/3){G{2L/3) - [G{L/3)} 2 } - [G{L/3)} 3 . Col- 
lecting terms, we obtain: Z2 = [G{L/3)] 3 {v — l) 2 + 
2[0(L/3)]0(2L/3)(v-l)+0(L). 

Let us define a "bubble" as a polymer segment with 
only its ends bound. The minimum length of a bub- 
ble is L m = L/{N +1) and the maximum length is 
equal to the total contour length of the filament L. No- 
tice that in our model a bubble will always have length 
< L <C L p . A "chain" is defined as a sequence of minimal 
bubbles (of length L m ). Fig. 1 provides a pictorial def- 
inition of "bubbles" and chains. The partition function 
Zn is a sum which consists of all products of the form 
f{m 1 )g{ni)f{m2)g{n 2 ).:f{m k ) where 



f(m) = 



V3 L P B 



(13) 



2tt L 2 n m 2 

is the statistical weight of a bubble of length mL m , 

3 (n) = [/(l)]>-ir +1 (14) 
is the statistical weight of a chain of length nL m , and 



mi+ni + m 2 + n 2 + ...+m,k = N+l, 0<k< 



N+l 



where [(TV 
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l)/2] is the integer part of {N 
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l)/2 . We 
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The first sum in Eq. ( Jig ) represents configurations with 
only one bond and we shall denote it by ■ The curly 
braces indicate that the sums must satisfy the constraint 
of Eq. ©. As N increases, calculating the combina- 
torial factors becomes an impossible task. That is why 
we use a standard trick and incorporate the constraint in 
the partition function via a Kronecker delta j2j| : 

[(JV+l)/2] oo oo fe 

z n = d n + EE^-EK-+ m ^ x 

k— 1 n— 1 m— 1 j—l 

A Tit i /(w,)p(n,) 

exp{/3[iV - + "*j)]> J J w V - ( 17 ) 

J = 1 9{ n k) 



z-plane 



i \ ^ X \ ^ X \ ^ X \ ^ X \ ^ X \ ^ X 

wnere 2^ n =i 2^ m =i = 2^m=i •■■ 2^n fc =i 2^ mi =i ■■• 2^ mfc =i 
and the auxiliary real parameter /3 has been introduced 
to make sure that the partition function converges in 
later steps of the calculation. 

We now introduce a complex representation of the Kro- 
necker delta which yields: 
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where 
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(f>j and are independent of j, 



with z = exp[- 
Since both 

n j= i = ($^')' c , and since we are interested in 

a very large number of pinning sites (thermodynamic 

limit, N — > oo), we approximate Ysk=i~ (^^) h ~ 
^^/(l — <E >1 3>) . Using analytic continuation, we transform 
the integral over 8 to a contour integral over the complex 
"fugacity" z where the contour encircles the origin z = 
in the counterclockwise direction once and we obtain: 



Z N = D N + — I dzz' 1 *' 1 * 
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FIG. 2: The complex z plane with poles at z — 0, z — zo and 
a branch cut along the positive real axis starting at z = 1. 
The contour C can be deformed to a contour around zo and 
a contour around the branch cut which closes at infinity. 



with L 2 (z) being Euler's dilogarithm function. p4[ The 
integrand in Eq. (^]J) has three singularities: a pole of 
order TV + 1 at z = 0; a simple pole at the solution z of 
equation 



V(zoMzo) = 1; 



(24) 



and a branch cut along the positive real axis starting 
at Z = 1 due to L 2 {z) . The contour C encircles only 
the singularity at the origin because of the assumptions 
that we had made in deriving Eq. (|2l]). That is, \z\ 
was chosen so that the series $ and \& converge and also 
< 1 . As shown in Fig. 2, the contour C can be 
deformed into a contour which encircles only zq clockwise 
and a loop which goes around the branch cut and closes 
at infinity counterclockwise. [p5| At sufficiently low tem- 
peratures, < zq < 1 . In the thermodynamic limit, -D/v 
vanishes and the partition function is determined by the 
pole at zq: N^ 1 In Zn ~ — lnzo • 

As the temperature increases, v decreases and it can 
be seen from Eq. (^J) that zq is shifted to the right. In 
the thermodynamic limit, Eq. ( |ll|) yields Q = — [(v — 
l)/z ](dzo/dv). This implies that the average fraction 
of intact bonds, which is a physically observable quan- 
tity, monotonically decreases.]^ The unbinding transi- 
tion occurs when zq — ► 1 . Therefore, the unbinding cri- 
terion is $(1)*(1) = 1 . Given that L 2 (l) = C(2) « 1.64, 
where £(z) is Riemann's zeta function, we obtain 



h^[exp(J/k B T c 



2 . 



(25) 



This is the main result of the article. The transition 
will be of second order because the derivative of L 2 (z) 
diverges logarithmically at z = 1 and therefore the aver- 
age fraction of intact bonds vanishes continuously at the 
critical temperature T c . Notice that in order to obtain 
this transition we need the thermodynamic limit where 
N — > oo and 10-i%N 2 B which is the probability of find- 
ing a pinning site within the binding region remains con- 
stant <C 1 . This implies that B — > which is consistent 
with a microscopic effective length that characterizes the 
bond. Of course, real systems will have a finite number 
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of pinning sites and the transition will not be sharp. For 
a sufficiently high density of pinning sites, however, we 
would expect a clear crossover from a low temperature 
phase with most of the pinning sites bound to a high tem- 
perature phase with most of the pinning sites unbound 
which will be described by Eq. (|25|). 



IV. BENDING RIGIDITY AND 
FORCE-INDUCED UNBINDING 

In this Section we consider the force- induced unbinding 
of a weakly bending semiflexiblc polymer and we show 
that the bending rigidity facilitates the unbinding. If we 
apply a transverse force to the free end of a clamped 
semiflexible polymer, the effective free energy of Eq. (Q) 
changes by an extra term: 



ds 



8t(a) 



ds 



/ / dssm9(s) 
Jo 



(26) 



where, as in Sec. II, 8(s) is the slope of the tangent vector 
with respect to the longitudinal direction. The partition 
function of this system is a path integral over all possi- 
ble conformations. Slicing the i-length into N segments 
each of length a and using the small-angle approximation, 
we obtain: 



N 



2—1 



fa 



N 

E< 

i=l 



(27) 



In this approximation, the path integral is Gaussian and 
can be easily calculated jtq] yielding 



Zf = exp 



iL p {k B Tf 



The corresponding free energy is 

f 2 L 3 



F = - 



3L p (fcsT) 



(28) 



(29) 



Notice that this free energy is just minus the elastic en- 
ergy of a cantilever spring, U e i = / 2 /2%, whith spring 
constant % — 3k/£ 3 . The latter has been obtained in a 
linear response calculation in Ref. [p7| . 

Although the original "Hamiltonian" (Eq. |2^) is ex- 
tensive, the free energy F of Eq. ( p9|) is not (grows as L 3 ) 
because Eq. ( |29| ) is an approximation valid only in the 
short length scales of the weakly bending limit (L < I,,). 
Given that the free energy of the bound state is always 
extensive, this non-extensivity leads to a "lever-arm ef- 
fect" in the force-induced unbinding. That is, for a long 
enough total length, the unbound state will be favorable 
having a lower free energy. We can estimate an upper 
bound for this "critical" length. If T\> is the free energy 
density of the bound state, the "lever-arm" critical length 
Li should satisfy the condition: 



f 2 L? 



3L p (k B T) 



(30) 



If L > Li, the transverse force / will always unbind the 
polymer (in equilibrium). Using the model of Sec. Ill, 
we have T h « {k B T/L m )[J/k B T - hx(V3L p B / L^tt)] . 
It turns out that 



Li 



*7 



J — fcnTln 



( V3L p B 
V LI 2ir 



1/2 ( ^ } i/2 

Lm. 



(31) 



For pulling forces of the order of picoNewton, persis- 
tence length of the order of /j,m, and binding free en- 
ergy per L m of the order of k B T, it turns out that 
Li « 10- 3 {L p /L m )^ 2 L p which is an indication of the 
relevance of the "lever-arm effect" to biopolymers. 

This is a phenomenon related to the "molecular lever- 
age" discussed in Ref. Q. In both cases the bending 
rigidity facilitates the force-induced unbinding. The two 
phenomena, however, are different. We describe a sit- 
uation where a bound state with a sequence of pinned 
sites becomes thermodynamically unstable when the sys- 
tem is long enough for the "lever-arm" to dominate the 
free energy whereas Ref. [^8| presents an estimate of the 
torque induced force exerted on a single ligand-receptor 
pair which turns out to be much stronger compared to 
that applied in traction. 



V. DISCUSSION AND CONCLUSIONS 

As we mentioned in the Introduction, the unbinding 
of semiflexible polymers is a particularly tricky problem 
because of the ambiguity of the relevant thermodynamic 
limit and also because of a lack of exact solutions of the 
WLC model with a binding potential. Refs. p, 113, ac- 
tually deal with the unbinding of flexible polymers(with 
L ^> L p ) and consider the effect of the bending rigidity 
on the conformational properties of the adsorbed (low- 
temperature) phase. In these works, the bending rigidity 
enters as a perturbation to the flexible (Gaussian) chain 
and the inextensibility is absent. The early references 
H % use a somewhat inconsistent approach where the 
scaling behavior of a weakly bending stiff filament (valid 
only for L <C L p ) is artificially extended to apply to any 
length. Another serious drawback of this approach is 
that it yields results which appear to be independent of 
the persistence length while, in principle, they should 
not. The idea is to solve Eq. (||) for large L with a 
binding potential using scaling Ansatze and then invoke 
the necklace model to predict the order of the unbinding 
transition from the scaling behavior of G. The details of 
the necklace model ( "bubbles" , "chains" , partition func- 
tion, etc.), however, are not worked out. Ref. J8| em- 
ploys a discrete model for stiff filaments where an extra 
("monomer") length is introduced and turns out to be 
relevant for the unbinding transition. It also proposes 
an energy/entropy melting criterion which applies to ad- 
sorbed phases similar to those discussed in |9], If . Ref. 
p2|, using a RG treatment, demonstrates the relevance 



of an orientation-dependent interaction field for the un- 
binding transition. The RG flow, however, implies a ther- 
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modynamic limit which carries on the inconsistencies of 
Refs. H 0. Ref. Jl3| models a semifiexible polymer as a 
directed self-avoiding random walk and it reiterates the 
inconsistencies of Refs. || |?J because the unbinding tran- 
sition occurs at the thermodynamic limit of an infinitely 
long walk (infinitely longer than the persistence length) 
where one would normally expect to recover the behavior 
of a flexible chain. 

In conclusion, applying the necklace model we have 
obtained a criterion for the depinning of anchored semi- 
flexible polymers in the weakly bending (stiff) limit. This 
model has been extensively used to study the unbinding 
of flexible polymers g |l| |j| This is its first de- 
tailed application to the unbinding of semifiexible poly- 
mers. A general and rigorous theoretical treatment of 
the unbinding of semifiexible polymers of arbitrary to- 



tal length and persistence length in the presence of an 
arbitrary binding potential has not yet been achieved. 
Our model suggests an alternative way to consider the 
thermodynamic limit for this system and straightens out 
several misconceptions of previous studies. 

We have also shown how the bending rigidity facilitates 
the force-induced unbinding of semifiexible polymers in 
the weakly bending limit. We have estimated a critical 
length as a function of the pulling force, the binding free 
energy density, the persistence length and the density of 
pinning sites above which the polymer acts as a lever-arm 
and unbinds. 

We thank T. Franosch, G. Gompper, M. C. Marchetti 
and A. Parmeggiani for useful comments and discussions. 
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